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Abstract — The main objective of the paper is to review 
a number of widely used multilateral index numbers for 
international comparisons of purchasing power parities 
(PPPs) and real incomes that can be derived using the 
stochastic approach. The paper discuss that price index 
numbers from commonly used methods like the Ikle, the 
Rao-weighted, and an additive multilateral system are all 
estimators of the parameters of the country-product- 
dummy (CPD) model. The paper also presents the 
method of moments (MOM) as an approach to estimate 
PPPs under the stochastic approach and shows how the 
Geary-Khamis system of multilateral index numbers is a 
method of moments estimator of the parameters of the 
CPD model. 
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multilateral index numbers - method of moments - country 
product dummy (CPD) model 
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I. Introduction 

The stochastic approach, used in derived the purchasing 
power parities (ppps) which are essentially spatial price index 
numbers that provide measures of price level differences 
across countries or regions within a country. 

The Country Product Dummy (CPD) model represented a 
simple regression approach to measures price level 

differences in different countries or regions within cl country. 
The CPD model proposed by Summers (1973), as a model for 
filling missing price data in the context of international 
comparisons. Subsequently the CPD model used in the first 
few phases of the International Comparison Project (ICP) 
mainly as a tool for aggregating price data under the basic 
heading level. 

Rao (1995), proposed Weighted- CPD as a link between the 
CPD model and the geometric variant of the Geary-Khamis 
model. In addition, it considered a generalization of the CPD 
model, which used for aggregation above the basic heading 
level. 

In p. .= 7i D + n D +...+ ji D + tj D + r\ D*+...+ n D + u.. (1) 
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M and N is number of countries and commodities 
respectively. 

I, v are random disturbance terms which are 

ij S 

independently and identically distributed. 


II. The Country Product Dummy (CPD) Method: 


Rao (2004), used the (CPD) model to represent multiple 
regression approach to measure price level differences in 

different countries. The CPD model postulates that the Pjj , 

is the product of three components: 

p. ■ = p. • PPP . • V.. (2) 

ij i J ij 

Where: 

pj is the commodity i' th price, PPP j ' s the purchasing 
power parity is the price level of the i" th commodity relative to 
other commodities price is p and V ’S are a random 

ij 

disturbance terms which are iid distributed. The parameters of 
the model (pppj and pi) can be estimated: 

In p.. = In p. + In ppp . + In V.. (3) 

y i j y 


lnp.. = ri. + 7i. + u.. 

y i j y 


( 4 ) 


Where, Tj and 21 _ are defined as before and U are 

i J ij 
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random disturbance which are iid normally distributed 

u : N ( 0 , a 2 ). 

ij 

It is possible to express Pj relative to a reference 
country (say country 1), then it j represents the purchasing 


power parity PPP j of country j which showing the number 

of country j currency units that have the same purchasing 
power as one unit of currency of country 1. The PPP for 
country, j is given by: 


PPPj = exp 




(5) 


This system solved by imposing a linear restriction 
on the unknown parameters. For a.i = 0 is the restriction 
imposed, it can be easily shown that, for each j=2,..., M 


pppj= exp(Tij) 


N 

n 

i=l 


Pij 


-|1/N 


( 6 ) 


Using Eq(6),to make a comparisons of price levels 
between two countries c and d, represented by ppp c( j can be 
derived as: 


PPPcd 


PPPd 

PPP n 


N 

= n 

i=l 


Pid 

p ic 


-il/N 


The price level comparison, PPP cd , in Eq(7) satisfied 
the transitivity property and Eq ( 7) is equivalent to the EKS 
(Elteto-Koves-Szulc) index number. 


The restriction ji = 0 imposed in Eq (8), this restriction 

implies that the currency of country 1 taken as the numerator 
and PPPs of all other countries are expressed relative to 
the currency of country 1. The model under this restriction is 
essentially the same as in Eq (9) except that the first column of 
X and the first element of P are dropped from the model. The 
restricted model is then given by: 


y-X*0*=v (10) 

Under the assumption of independently and iid random 

disturbances, the best linear unbiased estimator of (3 and 
the associated covariance matrix is given by: 


*r 


p = X X 


x*' y 


( 11 ) 


Var 



= q2 


X X 



( 12 ) 


Using the special structure of the matrix X*, consisting of 
various dummy variables as columns, it can be shown that: 


*' * 
X X 


(7) 


N Vl)x(M-l) V-1 )xN 

^Nx(M-l) M ^NxN 


(13) 


Where I is an identity matrix and J is a rectangular matrix of 
appropriate order whose elements are all equal to 1. Using the 
v*' v* 

special structure of A A : 


III. The CPD Regression Model 

The standard CPD formulation is multiple regression model 
that regresses logarithm of observed prices on a set of dummy 
variables representing the commodity and country to which a 
given price observation refers to as in the Eq (1) this model 
can be written as: 


y-- = 

X. . 

P + u 

00 

u 

y 

y 

Where: 


* * 

* 

x.. = 

U 

D i 

,D 2 ,„. ,D m , D t ,D 2 ,.. 

--’ d n 

P = 

7 T^ , 

71 Nr ’ yi 5 y 2 ’• 

— 


The values of the dummy variables are determined at the 
observation ij of NM observations ( for j =1,2,...M and i 
=1,2,...N). 


The model can be rewritten as: 


* 1 * 
X X 


i t i r l , 

N A (M-1)x(M-1) + N J (M-1)x(M-1) NM -(N+M) (M-1)*N 


1 


■ 


NM - (N + M) Nx ( M -') 

ven this expression, it can be shown that: 


1, M-lr 

M N+ MN J NxN 


1 ? 
71 . = - X 

J N i=i 


Pij 

Pi 


Gi 


(15) 


Rao (2004), also proved that 71 j is an unbiased estimator of 
71 . and its variance is given by: 

J 


2 2 

V| 71: | = 


(16) 


To estimate the variances of these parameters, the least 
squares residuals defined for each observation y as: 


y - xp = u (9) 

Eq (9) is a general regression model with NM observations 
and (N+M) explanatory variables. The matrix X has rank 
(N+M-l) as the first M dummy variables and the last N 
dummy variables sum to the same vector with elements equal 
to 1. 




(17) 


Since, ^ jj ’s are assumed to be iid with mean zero and 
variance a 2 . 

An unbiased estimator of a 2 is given by: 
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2 

CT = 


N M ~ 

Z Zeg 


(18) 


MN-(M+N-1) 

It is possible to derive the estimated variance of each element 

/v 

of P . Eq (18) is illustrated that the estimated standard errors 
of the it .’s are the same for all j=l,2,.. .M. Therefore, there 

is little that distinguishes the estimators of the logarithmic 
price level differences in terms of standard errors. This is 
mainly due to the fact that all commodities are priced in all the 

countries and disturbances in equation, ^ jj are i id. 


IV. Systems of price Index Numbers for 
International Comparisons Based on Stochastic 
Approach: 

These systems are a class of index numbers for international 
price comparisons, for example Rao system, Geary-Ghamis 
system (GK), Hajargasht-Rao, Elteto-Koves-Szulc (EKS), 
and IKLE system. Diewert (2005) has demonstrated that a 
number of commonly used formulae can be derived using the 
CPD model. Although, Rao (2005) established that the Rao 
method for computing PPPs is equivalent to the weighted 
CPD model. 

Rao, (2005) supposed that P— and C[— are the 

price and the quantity of the i-th commodity in the j-th 
country respectively, where j =1,2,...,M indexes the 
countries and i =1,2,...,N indexes the commodities. Also 

defined PPP j as purchasing power parity or the general 

price level in the j-th country relative to a numerator country , 
P j is the world average price for the i-th commodity and 

supposed that the following system of weights W.. and 

* 

w .i is defining different systems of index numbers. These 


weights are defined as: 


W 

ij N 


P ij q ij 


i - | l ’i.i q U 


W 


w 


1.) 


M 

X W 

j=l u 


Z w „ =1. 


5X=i. 


(19) 


( 20 ) 


Lemma(l): 


The Rao system is the weighed geometric averages in 
the place of arithmetic averages used in the GK system as 
follow: 


M 


\W 


PPP, 

V JJ 


N 

PPP;= n 


fr, A w ;: 


J 


j = l 


v P i y 


i =1,2,...,N 


j = l,2,...M 


( 21 ) 


( 22 ) 


These two equations of the Rao-system are equivalent to 
the weighted least squares estimate for the weighted CPD 
model. 

V. HAJARGASHT AND RAO MULTILATERAL SYSTEM OF INDEX 

Number: 

Hajargasht and Rao (2010), proposed their index number 
system , (H-R) with expenditure share weighted, it is 
arithmetic averages. The system is: 


P: 


PPPj = 


r M p ii 

Z — Wii 
j=lpppj J 


N p ij 

vi=l Pi J 


i=l,2,...,N 


j =1,2,...,M 


Where, , PPPj 


p.. , W.. , W 
1J 1J 


(23) 


(24) 


and D. 


1J 


Dj are defined as before. 


Lemma(2): If the disturbances terms LI — distributed 

gamma (r,r) in the derivation of Hajargasht and Rao 
multilateral system with expenditure share weighted 
arithmetic averages, then the H-R multilateral system is 
equivalent to the weighted maximum likelihood estimators of 
parameters of the CPD model. 


VI. IKLE Multilateral System of Index Number: 

The IKLE (1972) system, is an additive consistent 
system that is similar to the GK system and makes use of the 

two concepts of PPPj anc * international prices p - : 


PPPj 


p i 


N Pi 
Z Wj. 
i=lPii J 


M PPP; * 

Z - -W.. 

IJ 


j=l p 


IJ 


j=l,2,...,M 


i = 1,2,...,N 


Where: p. , ppp. , p.. , W.. , W 

*1 ^lj ij 


(25) 

(26) 
and D j 


IJ 


Dj* are defined as before. 
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Lemma(3): Hajargasht and Rao (2010), proved that if the 
U.. disturbance term distributed gamma (r,r) the IKLE 

1J 

system of index number is the MLE of the parameters from 
the CPD model take the form: 

1 1 u.. (27) 


>] 


PiPPPj 


ij 


Where U-j are random disturbance terms which are iid and 
as before they are assumed to follow a gamma distribution. 


VII. The Geary-Khamis Multilateral System of 
Index Number: 

The Geary-Khamis (Geary, 1958; Khamis, 1970) 
multilateral system, is a system of index number which based 
on the same two concepts PPP j °f currencies, and averages 

of prices - p - . Computations based on a simultaneous 

equation system as follow: 

N 

Z Pijqpj 

PPPj = (28) 


X P iqij 


i=l 


M / \ 

.ypijVpppj) 

P i = 11 -M- 


M 

X qjj 
j=l 


(29) 


Where: P- 


PPPj 


P • • and Q.. defined as 

IJ IJ 


before 


Haj argasht and Rao (2010), proved that the G-K 
multilateral system of index number is equivalent to the 
method of moments (MOM) estimator of the parameters of 
the CPD models as follow. 

* 


p.. = p. ppp. u.. 
ij i J ij 


(30) 


* 

Where U.. are random disturbance terms, which are iid, also 

IJ 

assumed that Eq (30) can be written in the following 
equivalent form. 


Pij 


Pi PPPj 


-1 


u.. 

ij 


(31) 


With E ( U.j ) = 1 the form of Eq (31) is a non-additive 

nonlinear regression model. 

A. Estimation of non-additive nonlinear models: 

Hajargasht and Rao (2010) considered that the CPD 
model is a non-additive model and then used the method of 
moment estimation technique to estimate the parameters. 
Assume the nonlinear regression model: 


r(y.,x. ,P) = u. , i=l,2,...,N (32) 

Where: 

y. is the dependent variable, U. is the random errors, 

1 

X. is a L x 1 vector and (3 is a 1 x K column 

1 

vector, r(y., x. ,(3 ) is a nonlinear function and i 

=1,...., N indexes the number of observations and also 
assumed that E( ) = 0 . The starting point to base the 
estimation of parameters in on the moment conditions is : 

E (X' u) = 0 (33) 

Where, X is NxL matrix containing X. S and u is an Nxl 
vector containing U j 8 

Generally the estimation based on the following K 
moment conditions: 

E(R(x,P)'u) = 0 (34) 

Where: R is a Nx K vector of functions of X and P . 


iy construction there are as many moment conditions as 
larameters. Therefore, a method of moment estimator can be 
>btained by solving following sample moment conditions: 

~(R(x,pyr(y, x,P)) = 0 (35) 

rhis estimator is asymptotically normal with variance matrix: 

(36) 


/here: 


i-l 


dr (y, x ,p ) 


IN , R = R ( X,p ) 


„ 2 U U 
a - 




B. Estimation ofPPPs under the optimal choice of moment 
conditions and standard errors using MOM: 

To postulate the observed price of j-th commodity in i-th 
country, pjj is the product of three components, the 

purchasing power parity PPP ,the price level of the j-th 

commodity relative to other commodities, p. and a random 

disturbance term as in Eq (30) as: 

p.. = p. ppp. u.. 

IJ 1 J IJ 

* 

Where: U.. are random disturbance terms, which are iid. 

IJ 

Rewrite this model in the following equivalent form as in Eq 
(31) as: 


P ij 


Pi PPPj 


-1 = u.. 
U 


With E ( U.j ) = 1 this equation is a non-additive 

nonlinear regression model. This model to be solved can be 
written as: 
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To impose weights in the price index, the R matrix will be 
(38) defined as: 


Where, R is an (M + N )( M + N ) matrix and it can 
be shown that most efficient choice of R is defined as: 


R'=E 


P„ 


PlM 


Pm 

PiPPPi 


PlM 

Jpppm 


Jim 

iPPPm 


Pnm 

»PPPm 


Pm 

nPPPm 


0 

Pi 

0 


PPP1 


q 12 . N2 

PPPi PPPi 
0 0 0 


_12_ q 1M 

Pi 

.% 

Pn 

0 0 


Wu 


W NM 


PPPm PPPm 


(42) 


Then the system of normal equations has the following 
forms: 


Since, 


ij 


Pi PPPj 


— 1 and f.. — 

“ 1 IJ 


IJ 


Pi PPPj 


-1, then: 


R'r: 


.1 

0 

.1 

0 

1 

0 


p» -1 


0 

Pi 


Pi 


Pn 



Pj PPPI 




1 


1 


1 




(1 
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0 



Pll 1 




Pn 


Pn 
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P2PPP1 



1 
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0 

0 
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PPPi 

PPPi 
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pppi 
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1 
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Pnm 1 



0 

0 
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PPPm 

PPPm 
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oj 


Then the system of normal equations has the following 
forms: 


i M 

Pi = - 2 
M J=1 


1 N 

PPPj=— s 
J N i=l 


r lJ 


Pi PPPj 
Pi 


IJ 

Pi PPPj 


(41) 


M 

Pi = 2 

j=l 

N 

PPPj = z 

J i=l 


Pi 


11 * 

w 


PPPj ij 


-W; 


Pi 


(43) 


This set of equations is the same as equations that 
defined the new system based on the expenditure share 
weighted arithmetic means to define pppj’s and pi’s. 

C. Estimation ofPPPs in Gk system under the optimal 
choice of moment conditions and standard errors using 
MOM in: 


To derive weighted Geary-Khamis (GK) price index, the 
quantity weighted price index R can be defined as: 


. jlM 

P, P N 

%i 

P, P» 


R’= 


PPP, 


.VL 

ppp, 


.Ml. 

PPP,' 


PPP, 


. Vl 
PPPm 


PPPm 


Then the system of normal equations has the following 
forms: 
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M 

Pi = I 
j=l 


Pij q ij/ Pp Pj 


N 

PPP;=I 
J i=l 


M 

I qj 

j=i ‘ 


Ml 

M 

I q ij 

j=l 


(45) 


Which is identical to the equations of the Geary-Khamis 
system given in Eq (28) and Eq(29). Thus it is clear that the 
G-K PPPs and Pj’s are the method of moments (weighted) 
estimators of the parameters of the CPD model. 
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